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We study in this paper the perturbations of the quintom dark energy model and the effects of 
quintom perturbations on the current observations. Quintom describes a scenario of dark energy 
where the equation of state gets across the cosmological constant boundary w = —1 during evolution. 
We present a new method to show that the conventional dark energy models based on single k-essence 
field and perfect fluid cannot act as quintom due to the singularities and classical instabilities of 
perturbations around w = —1. One needs to add extra degrees of freedom for successful quintom 
model buildings. There are no singularities or classical instabilities in perturbations of realistic 
quintom models and they are potentially distinguishable from the cosmological constant. Basing 
on the realistic quintom models in this paper we provide one way to include the perturbations for 
' dark energy models with parametrized equation of state across —1. Compare with those assuming 

, no dark energy perturbations, we find that the parameter space which allows the equation of state 

to get across —1 will be enlarged in general when including the perturbations. 
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INTRODUCTION 



In 1998 two groups Q have independently discovered the accelerating expansion of our current universe based 
on the analysis of Type la Supernovae (SN) observations of the redshift-distance relations. In the framework of 

■ Friedmann- Robertson- Walker (FRW) cosmology, the acceleration has been attributed to the mysterious source dubbed 
QQ \ dark energy. The simplest candidate for dark energy is a small positive cosmological constant, but it suffers from the 
"^\|" . difficulties associated with the fine tuning and the coincidence problem J3, Q . The most popular alternative to the 

■ cosmological constant is the model of rolling scalar field-quintessence ^ Q . In most cases the quintessence equation 
. of state (EOS) w changes slowly with time and can be well approximated with a constant w with il' > — 1 In the 

■ early probes of new physics, cosmologists have assumed a cosmological constant as the new component [ iL M M 11( 1 [Tl| 
\ and later fitted directly to the dynamical quintessence models 4, 12], or used a constant i^ 0, H, Il3l Il4| where w 

■ was restricted in the region oi w > —1. In Ref. 15] the author firstly extended the fitting of dark energy to include 
w < —1 and found some mild preferences. The author constructed a toy model of rolling scalar field with a negative 

Q ■ kinetic term and called it phantom TE^. The model of phantom has some theoretical problems 16] and there have 
^ I been many attempts towards resolving them llTl . 

c/^ . The accumulation of the observational data 1 1 SlIT^ [2(1 [21L [2312^ has opened a robust window for probing the more 

■ detailed behaviors of dark energy. There have been many studies in reconstructing the evolution of its energy density 
[2^ or equation of state 25, 26, 27j as a function of the redshift. Various parametrizations of w as well as dark energy 

. . - models have also been considered to fit directly to the observational data (e.g. j28l. l29l l30L l3lL 133. 133. 133. l35l l36l l37| ^ . 
rN \ Based on the fact that current observations cannot exclude dark energy models with the equation of state getting 
^ ■ across —1 during evolution with the redshift, we proposed a model dubbed quintom 35]. The model of quintom is a 
new scenario in the sense that the conventional quintessence or phantom models cannot realize the crossing of the 
cosmological boundary. Along this line the author in Ref. 38] has demonstrated that in the framework of general 
relativity the model of k-essence 39], where the scalar field of dark energy has non-canonical kinetic terms, cannot 
realize such a crossing behavior. A toy model with two rolling scalar fields which have opposite kinetic energy terms 
can easily realize the transition and it can be regarded as the simplest quintom model 35, 40]. Recently a single 
field quintom model was proposed in 41] by adding higher derivative operators in the Lagrangian. In the simplest 
case such a model is equivalent to the two-field case as proposed in In addition, the quintom model of dark 

energy are different from the quintessence or the phantom in the determination of the evolution and the fate of the 
universe. Due to its distinctive properties, the quintom model with oscillating equation of state across —1 can lead 
to the oscillations of the Hubble constant and a new scenairo of recurrent universe 4^], which to some extent unifies 
the early inflation and the current acceleration of the universe. Recently there have been a lot of interests in the 
phenomenological studies relevant to quintom models in the literature 0? QE QE EE 1^ EE EE Ell • 

Current supernovae data alone, which make the only direct detection of dark energy, seem to favor a quintom-like 
model at around 2-cr level |^ EE E3 • The quintom model is also mildly favored in the combined analysis with the 
cosmic microwave background (CMB), large scale structure (LSS) and supernovae data "3^,^34^. However, when some 
other observational data sets (such as the new observational data based on Chandra measurements of the X-ray gas 
mass fraction in 26 X-ray luminous galaxy clusters or the recent new constraints from the bias and Lya forest 
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of the Sloan Digital Sky Survey (SDSS)) have been taken into account the situation changes and the preference for 
quintom-like dark energy models becomes weak [s^ls^. However the previous fittings in the literature on quintom-like 
dark energy models have either fully or partially neglected the perturbations, which in some sense do not describe the 
realistic models with EOS across —1 and will lead to some bias in the fittings. The aim of this paper is to develop 
a self-consistent way to include the perturbations of quintom in light of the observations. We will present a simple 
new method to show that conventional single perfect fluid and k-essence dark energy models cannot act as quintom, 
which is due to the singularities and classical instabilities of perturbations. Based on the realistic quintom models in 
this paper we will provide one way to include the perturbations for dark energy models with parametrized equation 
of state across —1. Compared with those assuming no dark energy perturbations, we find that when including the 
perturbations the parameter space which allows the equation of state to get across —1 will be enlarged in general. 

This paper is organized as follows: in section II we discuss the difficulty of quintom model buildings and provide a 
new proof regarding the impossibility of single perfect fluid and k-essence model as quintom, then present some viable 
quintom models; in section III we study in detail the perturbations of the quintom models; in section IV we investigate 
the possible signatures of quintom models of scalar fields and the effects of quintom perturbations on the observations; 
in section V we provide one way to include the perturbations for models of dark energy with a parametrized equation 
of state across —1; we conclude in section VI. 



II. QUINTOM MODEL BUILDING 

A. Difficulties of Quintom Model Buildings 

We start with a brief overview on the arguments against the possibility of realizing the quintom scenario with a 
single fluid or a single scalar field in the conventional framework. 

Consider firstly a single perfect fluid, the energy-momentum tensor has the conventional form, 

T^z. = -P9^,v + (p + P)u^u^ , (1) 

where p and P are proper energy density and pressure, is 4- velocity with u^u^ = 1. The energy density and the 
pressure of the fluid can be parametrized as J5] 

P = /H , 

P = nf{n) - f{n) , (2) 

where f{n) is a positive function of n. The introduced variable n can be identified with the number density and 
the prime represents derivative with n. The equations of motion are just the covariant conservation equations of the 
momentum tensor, V ^T^^ = 0. In spatially flat FRW spacetime 

ds^ = a^{r){dr'^ - dx'dxi) , (3) 

there is only one equation, 

p + 37Y(p + P) = , (4) 

where the dot is the derivative with the conformal time r and Ti = d/a. Combining Eq.(4) with Eq. Q, we get 

/(n)(n + 3Hn) = . (5) 

Since f'{n) does not vanish everywhere (otherwise it corresponds to the cosmological constant), one has the conser- 
vation equation of the number density h + SHn = 0. In the expanding universe, n will decrease monotonically with 
time. 

In the following we will demonstrate that the system suffers from the problem of singularity and classical instability 
when the equation of state of the perfect fluid crosses the boundary of —1. Let us assume that the system crosses 
— 1 at the point of n = no 7^ 0. At this point p{no) + P{no) = 0, /'(no) = and /'(n) will change the sign after the 
crossing. So, in the neighborhood of no, we can expand f'{n) in terms of (n — no). The adiabatic sound speed square 
in this neighborhood is 

dp f'{n) ^ n-no' ^ ' 
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We can see that is singular at the crossing point. Moreover, is negative in the region of n < uq. And when n 
approaches no from this side, it wih approach — oo. A negative wih induce a serious classical instability to the system, 
the perturbations on small scales will increase quickly with time and the late time history of the structure formations 
will get significantly modified. This will inevitably lead to the fact that such models can never be compatible with 
the observations relevant to structure formations, such as CMB and LSS. 

As shown above it is impossible to realize the quintom scenario with a single perfect fluid, now we turn to consider 
the model of scalar field. The general model of dark energy with a single scalar field and an arbitrary function of its 
first derivative in the Lagrangian was proposed in Ref. 39] and named as k-essence. Its Lagrangian usually has a 
non-canonical form 

C = P{cP,X) , (7) 

with 

X = ^V^<PW4> ■ (8) 
The energy- momentum tensor of this system has the same form as that of the single perfect fiuid Eq. where 



P = 2XP,. 



(9) 



Let's see under what conditions the system will be able to cross the barrier of w = —1. In order to do that, one 
requires p + P to vanish at a point of (^o, -^o) and change the sign after the crossing. This can only be achieved by 
requiring P,x (^o, -^o) = and P,x has different signs before and after the crossing since X cannot be negative. The 
covariant conservation law of the energy- momentum tensor gives the equation of motion, 

{P,x 5^" + P,xx W<PV^<P)V^W,<P + p,^=0. (10) 
From this equation, we obtain the equation for the background field 

p,x [4> + (34 - IW] + = , (11) 
and the perturbation to the first order (neglecting the metric perturbations for the time being): 

u + [-4V2 - z/z - 34(W - H')]u = , (12) 

where 



h2 



and the effective sound speed is given by 



P.x I , (13) 



(14) 



This sound speed c^^ is often used in describing the perturbations of the scalar fields instead of the isentropic sound 
speed, which behave differently due to the intrinsic properties of the scalar fields 56]. For a conventional quintessence 
or phantom field, c^^ = 1. The dispersion relation from Eq. (|T^ is 



0? = cl,k^ - z/z - 3cl,{n - U^) . (15) 

One of the conditions for the stability of k-essence perturbations is that c^^ must be positive ^3^ . This requires that 
p,x has the same behavior as that of P,j\:, i-e., it must vanish at the crossing point and change the sign after the 
crossing. 

Similar to the analysis in the case of single fluid, we can see that z/z diverges at the point (^o, -^o)- This singularity 
is unavoidable in the perturbation equation and the physical quantities describing the fluctuations are not well deflned. 
Generally, z does not vanish at the crossing point, hence there exists a region in which o;^ < and the perturbation 
is unstable. Furthermore, the canonical momentum deflned by the Lagrangian {TJ is 

U='^ = P,A. (16) 
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Its derivative with respect to ^, 



(17) 



vanishes at the point of crossing. This shows that (j) is not a single valued function of the momentum IT and we cannot 
get a canonical Hamiltonian transformed from the Lagrangian unambiguously The theory cannot be quantized 
in a canonical way. Hence we have shown that the conventional k-essence model cannot give rise to w across —1. A 
different proof is given in Ref. [23| . 

We should stress again that in realistic quintom model buildings one must consider the aspects of perturbations, 
where there are often dangerous instabilities in the conventional case. The concordance cosmology is based on the 
precise observations where many of them are tightly connected to the growth of perturbations and we must ensure 
the stability of perturbations. If we start with parametrizations of the scale factor 58j or EOS to construct quintom 
models, it can be realized arbitrarily if we do not consider the stability of perturbations. On the other hand when 
we start from scalar fields and use some phenomenological parametrizations it is in some sense very easy to resemble 
fluid behavior in the background evolutions. However the stability of perturbations must be considered. 



B. Some Viable Quintom Models 

As we demonstrated above in the conventional cases with a single fluid or a k-essence one cannot realize a viable 
model of quintom, we need to introduce extra degrees of freedom to realize the transition of w across —1. One of the 
possibilities is a system including two fluids with one being w > —1 and another w < —1. Speciflcally, consider a 
model which consists of two Chaplygin gases 59] with Pi = —Xi/pi and P2 = — A2/P2, in which Ai and A2 are positive 
constants. If pf > Ai and P2 < A2, one has > k;i > — 1 and W2 < —1. This system will cross the boundary of —1 
at some time because p2 is always increasing and pi decreasing. The sound speed squares are c^^ = —Wi > with 
i = 1, 2, hence the system will be free of the difliculties associated with the singularity and the classical instability 
which exist in the model of a single fluid. Furthermore, the flnal state of this system will be characterized by = — 1, 
the universe will approach the de Sitter space in the far future. In such a scenario there will be no big rip. In the 
framework of the fleld theory, the simple way to introduce the extra degree of freedom for the quintom model is the 
double scalar flelds model with one being quintessence-like and one phantom-like. We should point out that when 
adding extra degrees of freedom in the above way, this does not help solve the cosmological constant problem and 
nor can it help solve the coincidence problem, since for the component where < — 1 it cannot have the property of 
tracking behavior and has to be flne tuned ^. The above way of introducing more components provides the simplest 
possibility of quintom model building. 

There is another possibility of introducing the extra degrees of freedom for the realization of the transition from the 
quintessence phase to the phantom phase. This is the model proposed in Ref. 41] by introducing higher derivative 
operators to the Lagrangian. Speciflcally in 0| we considered a model with the Lagrangian 



>C = -^V^</'V''<^+^n<^n</'-n'/') , (18) 



where □ = V^V^ is the d'Alembertian operator. The term related to the d'Alembertian operator is absent in the 
quintessence, phantom and the k-essence model, which is the key to make the model possible for w to cross over —1. 
We have proven in this Lagrangian is equivalent to an effective two-fleld model 



with the following deflnition 



1 1 

£ = - 2 V^^/'V^V' + 2 V^XV^X - Vi^P - x) - ^X' , (19) 

X = -^^4> , (20) 

^ = 'P + X- (21) 

Note that the redefined fields xp and x have opposite signs in their kinetic terms. One might be able to derive the 
higher derivative terms in the effective Lagrangian (18) from fundamental theories. In fact it has been shown in the 



For the k-essence field where w < —1, it also needs to be fine tuned. 
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literature that this t ype o f operators does appear as some quantum corrections or due to the non-local physics in the 
string theory ,6p|. [6ll l6^ . With the higher derivative terms to the Einstein gravity, the theory is shown to become 
renormalizable which has attracted many attentions. In fact the canonical form for the higher derivative theory 
has been put forward by Ostrogradski about one and a half century ago 64]. In short, it is interesting and worthwhile 
to study further the implications of models with higher derivatives in cosmology (for a recent study see e.g. ^52j). 



III. PERTURBATIONS OF THE QUINTOM MODEL 

The quintom scenario as we have argued above needs extra degrees of freedom to the conventional models of a single 
scalar field, such as quintessence, phantom and k-essence and the simple realization of the quintom is a model with 
two scalar fields or "equivalent ly" two scalar fields for the case with the higher derivative operators. In the discussions 
below on the perturbations we will restrict ourselves to the two-field model of quintom with the following lagrangian: 

C = Cq^Cp (22) 

where 

CQ = \d^(^id>'(^i-V,{4>i) (23) 

describes the quintessence component, and 

Cp = -^d^hd^h - V2{h) (24) 

for the phantom component. The background equations of motion for the two scalar fields 0i(i = 1, 2) are 

dV 

^,^2n^,±a^--^ = , (25) 

where the positive sign is for the quintessence and the minus sign for the phantom. In general there will be couplings 
between the two scalar fields, here for simplicity we neglect them. 

For a complete study on the perturbations, the fluctuations of the fields as well as those of the metric need to be 
considered. In the conformal Newtonian gauge the perturbed metric is given by 

ds^ = a^{r)[{1^2^)dr^ -{l-2^)dx'dxi] , (26) 

Using the notations of ^651, the perturbation equations satisfied by each of the components of the quintom model (22) 
are: 

SP 

5i = -{l^w,){0,-3^)-3n{-^-w,)5i , (27) 

opi 

0, = -n{l-3w,)0,--^0,^k\^-p^S,-a,^^) , (28) 



where 



(^V<?^^)¥^, ^z=0 , (29) 

(30) 



pi 



and 



ucpi — upi -r 

Combining Eqs. ((^Tj), (|2H|) and we have 



SPi = 8 Pi - 2Vl8(^, = Sp, + ^[3H(1 - w^) + m] . (31) 



s^ = -(1 + w,){o, - 3^) - 3n{i - w,)6i - 3n^^^^^^^^—^e, , (32) 
Oi = 2ne,^- s.^kH. (33) 
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Since the quintom model in (22) is essentially a combination of a quintessence and a phantom field, one obtains the 
perturbation equations of quintom by combining the equations above together. The corresponding variables for the 
quintom system are 



Ei p 

and 

n T^iiPi^PiW^' 



'-'quintom — 9 v^^/ 



^quintom 



EM + Pi) 



(36) 



Note that for the quintessence component, —1 < wi < 1, while for the phantom component, W2 < —1. 

With the two fields the quintom model in (22) will be characterized by the potential Vi in (23) and (24). In this 
paper we take Vi{(j)i) = ^m|(/)|. In general the perturbations of 0^ today stem from two origins, the adiabatic and the 

isocurvature modes. If we use the gauge invariant variable d = — instead of Si, and the relation <I> = ^ in 

the universe without anisotropic stress, the equations (p?^ and can be rewritten as, 

Ci = -| - aid + $ + g^i) , (37) 

di = 2Wei + A:2(3Ci + 4$) , (38) 

where 

111 ■ 

Ci = + mil - w,) = ao[ln(a6|p, +pi|)] . (39) 

i -|- Wi 

(a is the curvature perturbation on the uniform-density hypersurfaces for the a-component in the universe |6^ . 
Usually, the isocurvature perturbations of are characterized by the differences between the curvature perturbation 
of the uniform-0i-density hypersurfaces and that of the uniform-radiation-density hypersurfaces, 

S^r = 3(0 - Cr) , (40) 

where the subscript r represents radiation. In this paper, we assume there is no matter isocurvature perturbations, 
so = C^. Eliminating d in equations (pTf)) and (pTHj) . we obtain a second order equation for 6>i, 

0, + (C, - 2n)0, + {QH - 2H + e)0, = k\4^ + Q$) . (41) 

This is an inhomogeneous differential equation, the general solution to it is the sum of a general solution to its 
homogeneous part and a special integration. In the following, we will show that the special integration corresponds to 
the adiabatic perturbation. Before the era of dark energy domination, the universe is dominated by some background 
fluids, for instance, the radiation or the matter. The perturbation equations of the background fluid are, 

C/ = -Of/s , 

Of = -n{l - 3wf)0f + k'^[3wfCf + (1 + 3^/)^] . (42) 

From the Poisson equation 

- = ^ E ""(1 + + * + p^") ^ + + * + ^(^f) > (43) 

a 

we have approximately on large scales, 

$ -C/ - ^Of . (44) 
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Combining these equations above with H. = 2/[(l + 3ti)/)r], we get (note numerically Of ~ 0(fc 

Cf = -T^TT— — = Const. , 

So, we can see from Eq. that there is a special solution to it which is given approximately on large scales by , 

ef = Of , (46) 

and from Eq. (|38j) we have, 

Cf = Cf ■ (47) 

This shows that the special integration to Eq. has the meaning that it corresponds to the adiabatic perturbation. 
Hence, for the sake of isocurvature perturbations of 0^, we can only consider the solution to the homogeneous part of 
Eq. (gU, 

Oi + {d - 2n)0i + {dn -2H^ k^)Oi = . (48) 
These solutions are represented by 6>-*^ and Q^^. The relation between them is 

niso oov/jiso 

cr° = '-^^^ ■ (49) 

Since the general solution of d is 

Ci = cr'+cr = Cr+cr , (50) 

the isocurvature perturbations are simply 5*^^ = 3Q^^. 

In order to solve Eq. we need to know the forms of Ci and H as functions of time r. For this purpose, we 

solve the background equations ((^ . In radiation dominated period, a = Ar , H = l/r and we have 



/ -A 9 . AT / -A 
-mir ) + A2J_i/4(- 

and 



'^'[AiJi/,{-m,r') + J_i/4(-mir2)] , (51) 



h = T-^'\Mhi^{\m^T^) + i2/-i/4(|m2r2)] , (52) 

respectively, where A, and Ai are constants, Jjy(x) is the vXhi order of Bessel function and Iv{x) is the z^th order of 
modified Bessel function. Usually the masses are small in comparison with the expansion rate in the early universe 
mi <C 7Y/a, we can approximate the (modified) Bessel functions as Jv{x) ~ (c\ + C2X^) and Iv{x) ~ (c\ + C2X^\ 
We note that J_i/4 and /_i/4 are divergent when x ^ 0. Given these arguments one can see that this requires large 
initial values of and if A2 and A2 are not vanished. If we choose small initial values, which is the natural choice if 
the dark energy fields are assumed to survive after infiation, only A\ and A\ modes exist, so will be proportional to 
in the leading order. Thus, the parameters Ci in equation ((S^ will be Ci = 10/r (we have used \pi^Vi\ = ^i/<^^)- 
So, we get the solution to Eq. (gSJ, 

^so ^ ^-4 [2).^ cos(A:r) + A2 sin(A:r)] . (53) 

$1^^ oscillates with an amplitude damping with the expansion of the universe. The isocurvature perturbations Q^^ 
decrease rapidly. If we choose large initial values for and A2 and A2 modes are present, <j)i will be proportional 
to in the leading order and Ci = 0. Now the solution to Eq. is 

l9f ^ = r[Ai cos(A:r) + A2 sin(A:r)] . (54) 

01^^ will oscillate with a increasing amplitude, so Q^^ remains constant on large scales. 
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Similarly, in matter dominated era, a = Br'^ , H = 2/r, the solutions for the fields are 

B B 

01 = r~^[Bi sin(— miT^) + B2 cos( — mir^)] , (55) 
o o 

and 

~ B ~ B 

02 = r-^[Bi sinh(— msr^) + B2 cosh( — msr^)] , (56) 
o o 

respectively. We get the same conclusions as those reached by the above analysis for the radiation dominated era. If 
we choose small initial values at the beginning of the matter domination, we will get the isocurvature perturbations 
in decrease with time. On the contrary for large initial values, the isocurvature perturbations remain constant on 
large scales. This conclusion is expectable. In the case of large initial velocity, the energy density in the scalar field 
is dominated by the kinetic term and it behaves like the fiuid with w = 1. The isocurvature perturbation in such a 
fiuid remains constant on large scales. In the opposite case, however, the energy density in the scalar field will be 
dominated by the potential energy due to the slow rolling. It behaves like a cosmological constant, and there is only 
tiny isocurvature perturbation in it. 

We have seen that the isocurvature perturbations in quintessence-like or phantom-like field with quadratical po- 
tential decrease or remain constant on large scales depending on the initial velocities. In this sense the isocurvature 
perturbations are stable on large scales. The amplitude of these perturbations will be proportional to the value 
of Hubble rate evaluated during the period of infiation Hinf if their quantum origins are from infiation. For a 
large Hi^f isocurvature dark energy perturbations may be non-negligible and will contribute to the observed CMB 
anisotropy'6?,'6^. In the cases discussed here, however, these isocurvature perturbations are negligible. Firstly, large 
initial velocities are not possible if these fields survive after infiation as mentioned above. Secondly, even though the 
initial velocities are large at the beginning of the radiation domination, these velocities will be reduced to a small 
value due to the small masses and the damping effect of Hubble expansion. In general the contributions of dark energy 
isocurvature perturbations are not very large 69j and here for simplicity we assume Hinf is small enough that the 
isocurvature contributions are neglie; ible^ Thus we will concentrate on in next sections the effects of the adiabatic 
perturbations of the quintom model with two scalars considered in this paper. 



IV. SIGNATURES OF QUINTOM AND THE EFFECTS OF PERTURBATIONS ON OBSERVATIONS 

Based on the perturbation equations and we modify the code of CAME and will study preliminarily 
in this chapter the observational signatures of quintom. Throughout this paper we assume a fiat universe. In showing 
the illustrative effects for quintom we have assumed the fiducial background parameters to be = 0.042, = 
0.231, r^D^; = 0.727, where DE denotes dark energy and today's Hubble constant is fixed at Hq = 69.255 km/s 
Mpc~^. We will calculate the effects of perturbed quintom on CMB and LSS. 

In the quintom model we focus on there are two parameters: one is the quintessence mass and the other being the 
phantom mass. When the mass of quintessence is heavier than Hubble parameter the field will start to oscillate and 
consequently one will get an oscillating quintom. In the numerical discussions we will fix the mass of the phantom 
field to be mp ~ 2.0 X lO^^^m^^. We vary the quintessence mass with the typical values being tuq = 10~^^ rripi 
and 4 x 10~^^mpi respectively. We plot in Fig. 1 the equations of state as function of the scale factor for the above 
two sets of the parameters and their corresponding effects on the observations. One can see the obvious oscillating 
feature of quintom as the mass of quintessence component goes heavier. After touching the w = —1 pivot for several 
times, w crosses —1 consequently where the phantom part dominates dark energy. The quintom field modifies the 
metric perturbations: Sgoo = 2o?^^5gii = 2a'^^Sij and consequently contribute to the late time Integrated Sachs- 
Wolfe (ISW) effect. The ISW effect is an integrant of i> + ^ over conformal time and wave number k. The above 
two quintom models yield quite different evolving <I> + ^ as shown in the right panel of Fig^ where the scale is 
k ~ 10~^ Mpc~^. We can see the late time ISW effects differ significantly when dark energy perturbations are taken 
into account (solid lines) or not (dashed lines). 

ISW effects take an important part on large angular scales of CMB and on the matter power spectrum of LSS. For 
a constant EOS of phantom Ref.[7l| has shown that the low multipoles of CMB will get significantly enhanced when 



^ We assume in the next section when the mass of quintessence is larger by an order and oscillates during late time evolutions, the 
adiabatic condition still satisfies well. 
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dark energy perturbations are neglected. On the other hand for a matter hke scalar field where the equation of state 
is around zero, perturbations will also play an important role on the large scales of CMB, as shown in ref.^^l- Our 
results on CMB and LSS refiect the two combined effects of phantom and oscillating quintessence. Note that in the 
early studies of quintessence effects on CMB, one usually considers a constant Weff instead: 

f daQ(a)w(a) , , 

/ dan{a) ' ^ ^ 

however this is not enough for the study of effects on SN, nor for CMB when the EOS of dark energy has a very large 
variation with redshift, such as the model of oscillating quintom considered in this paper. 

To face the oscillating model of quintom with the current observations, we make a preliminary fitting to the 
first year Wilkinson Microwave Anisotropy Probe (WMAP) TT and the TE temperature-polarization cross-power 
spectrum as well as the recently released 157 "Gold" SN data|23|. Following Refs.fT^. in all the fittings below 
we will fix r = 0.17, Qrnh'^ = 0.135 and Q^h'^ = 0.022, we set the spectral index as ns = 0.95 and the amplitude 
of the primordial spectrum will be used as a continuous parameter. In the fittings of oscillating quintom we've fixed 
the mass of phantom to be mp ~ 6.2 x 10~^^mp^ Fig |21 delineates 3a WMAP and SN constraints on the two-field 
quintom model, it also shows the corresponding best fit values. In the labels tuq and mp stand for the mass of 
quintessence and phantom respectively. The left panel of Fig|21 shows the separate WMAP and SN constraints. The 
green(shaded) area is WMAP constraints on models where dark energy perturbations have been included and the blue 
area(contour with solid lines) is without dark energy perturbations. The perturbations of dark energy have no effects 
on the geometric constraint of SN. The right panel shows the combined WMAP and SN constraints on the two-field 
quintom model with perturbations (green/shaded region) and without perturbations (red region/contour with solid 
lines). We find the confidence regions do show a large difference when the perturbations of dark energy have been 
taken into account or not. 

So far we have investigated the imprints of oscillating quintom on CMB and LSS. Now we consider another example 
where w crosses —1 smoothly without oscillation. It is interesting to study the effects of this type of quintom model 
with its effective equation of state defined in (57) exactly equal to —1 on CMB and matter power spectrum. This 
study will help to distinguish the quintom model from the cosmological constant. We have realized such a model of 
quintom in the lower right panel of FigO which can be easily given in the two-field model with lighter quintessence 
mass. In this example we have set tuq ~ 2.6 x 10~^^mpi,mp ^ 6.2 x 10~^^mp^ We assume there are no initial 
kinetic energy. The initial values of the quintessence component is set as = 0.226mpi and the phantom part: 
02i = 6.64 X 10~^ rripi. We find the EOS of quintom crosses —1 at z ~ 0.15, which is consistent with the latest SN 
results. 

The model of quintom, which is mainly favored by current SN only, needs to be confronted with other observations 
in the framework of concordance cosmology. SN making the only direct detection of dark energy, this model is most 
promising to be distinguished from the cosmological constant and other dynamical dark energy models which do not 
get across —1 by future SN projects on the low redshift (for illustrations see e.g. 32j). This is also the case for the 
model of quintom in the full parameter space: it can be most directly tested in low redshift Type la supernova surveys. 
In the upper left panel of Fig|31 we delineate the different ISW effects among the cosmological constant (red/light 
solid), the quintom model which gives Weff = — 1 with (blue/dark solid) and without (blue dashed) perturbations. 
Similar to the previous oscillating case, the difference is very large when switching off quintom perturbations and much 
smaller when including the perturbations. In the upper right panel we find the quintom model cannot be distinguished 
from a cosmological constant in light of WMAP. The two models almost give exactly the same results in CMB TT 
and TE power spectra when including the perturbations. We find the difference in CMB is hardly distinguishable 
even by cosmic variance. 

Given the fact aobve that from CMB observations quintom with Weff = — 1 makes no distinctive signatures, now 
we discuss briefiy the signatures in some other observations. To do that we need to consider the physical observables 
which can be affected by the evolving w sensitively. In comparison with the cosmological constant such a quintom 
model gives a different evolution of expansion history of universe, such as altering the epoch of matter-radiation 
equality. The Hubble expansion rate H is : 

H=^= Hoinma-^ + Qra-^ + X] (58) 

where X, the energy density ratio of dark energy between the early epochs and today, is quite different for the 
quintom-CDM and ACDM. In the ACDM scenario, X is simply a constant while in general for dark energy models 
with varying energy density or EOS, 

X = l]z,£;a-^e-'/"^^"^^^"" . (59) 
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The two models will give diflFerent Hubble expansion rates. This is also the case between the quintom model with 
Weff = — 1 in the left panel of Figl^land a cosmological constant. Different H leads directly to different behaviors of 
the growth factor. In the linear perturbation theory all Fourier modes of the matter density perturbations grow at 
the same rate. The matter density perturbations are independent of k: 

h + - AirGa'^pMSk = . (60) 

The growth factor Di{a) characterizes the growth of the matter density perturbations: Di{a) = Sk{a)/Sk{a = 1) 
and is normalized to unity today. In the matter-dominated epoch we have Di{a) = a. Analytically Di{a) is often 
approximated by the Meszaros equation g7^] : 

where we can easily see the difference between the model of quintom and a cosmological constant due to the different 
Hubble expansion rates. More strictly one needs to solve EQ. ((Sn|) numerically. In the lower left panel of Fig|31we 
show the difference of Di{a) between the quintom with Weff = — 1 and the cosmological constant. The difference in 
the linear growth function is considerably large in the late time evolution and possibly distinguishable in future LSS 
surveys and in weak gravitational lensing (WGL) observations. WGL has emerged with a direct mapping of cosmic 
structures and it has been recently shown that the method of cosmic magnification tomography can be extremely 
efficient 76j, which leaves a promising future for breaking the degeneracy between quintom and a cosmological constant. 



V. PERTURBATIONS OF PARAMETRIZED QUINTOM AND THE EFFECTS ON THE 

OBSERVATIONS 

There have been many studies in the literature in the fittings of the dark energy with parametrized EOS, such as 
the linear parametrization w = wq -\- wiz [T^ to SN and other observations such as CMB and LSS. For the latter 
observations the perturbations of dark energy need to be considered. However, at the point of w = —1, as pointed 
out in Section II one would be encountered with the singularity of the isentropic sound speed. Moreover in the 
perturbation equation one will get infinite 0. For the physical quantity {p -\- P)0 in the model of the single 
field of quintessence, it is not divergent at = — 1, ie. 6 ^ oo but {p -\- P)0 = k'^^dcj) = 0, however for the model 
with parametrized EOS one will generically have an unphysical divergence when it; 7^ at the cosmological constant 
boundary. The detailed explanation is given as follows: firstly from Eq.(|2HI) one will get infinite and the physical 
continuity implies that one will also get ^ 00 dit w = —1. Introducing the new physical quantity which is relevant 
to the CMB observations: 

V = {l^w)0 , (62) 

Eqs.(|S21 ESI) can be rewritten now as: 

^ = -V + (1 + w)3^ - mil - w)8 - 3wMi±^f!)±Mi^ V , (63) 
V = 2HV + i^^J + — ^V + /c2(l + ^)^ . (64) 

We can easily see that V ^ 00 when w ^ {) w = —1. 

We should point out that both the scalar fields and fluids obey the same form of equations on the evolution of 
perturbations: Eqs. (|27I28|) . and the only difference comes from the term of 5Pi/5pi. If one starts from Eqs.(|S21IS3 
and study the effects of dark energy by parametrizing the EOS, this is equivalent to the description of the effects of 
the scalar field and is identical to work starting with dark energy potentials. If in models with the parametrized EOS 
we have w always in the range [—1, 1], 01 w < —1 for < a < 00, there will be no unphysical divergence and this 
equivalently describes the single field of quintessence or phantom 78]. For example in model with w = wo-\-wi sin (In a), 
if we restrict wq = and \wi \ < 1 then Eqs. ((S51IM|) will always be continuous. However when the parameter space is 
enlarged to include w dit w = —1 Eqs. ((S51 IM)) will be unphysical. 

We emphasize that the above discussions are valid only for models with a single field. For models with multi fields 
we have shown explicitly in the previous sections the perturbation equations Eqs. ((S21IS5|) are continuous during the 
crossing of the cosmological constant boundary. It is similar for models with two fluids or models with two components 
of parametrized EOS: w = T^^liWi where each component Wi does not evolve across —1. This implies, however, in the 
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fitting of tiie models to tiie observational data the parameters introduced for the EOS should be doubled if allowing 
the EOS w to vary and get across —1. Certainly this is not practically applicable. It would be nice to develop a 
technique to include the perturbations which approximates well to the quintom, meanwhile not introducing the extra 
degrees of freedom to the models considered widely in the literature with parametrized EOS. We will make a proposal 
for it below. 

First of all we consider a system of quintom with two fields as above, being quintessence-like and ^2 being 
phantom-like, but restrict the EOS of the system not to cross over —1. In this case we will show the background of 
this system is equivalent to a model with an effective single scalar field denoted by x- By definition the pressure P 
and energy density p of the x field should be equal to the two-field case. When the kinetic term of 0i is larger than 
that of the phantom part the whole system of dark energy gives rise to an EOS larger than —1 and the effective 

• 2 -2 

X behaves like a quintessence. On the contrary when — 02 ^ 0, x is a phantom field. Hence the kinetic and 
potential terms of x, in terms of 0i and (/)2, can be expressed as 

± = <^i^ - <i>2^ (65) 

and 



V{x) = V{<P,) + V{<P2) , (66) 

where the "-h" sign in Eq. is for the case where the total EOS of dark energy is quintessence-like and the "-" 
sign for phantom-like evolutions. We can directly reconstruct the potential and time evolutions of x- For example if 
we set the potentials of the two fields to be both linear: 

^z(</>z) =^Oz + A,(/), , (67) 

in the early epochs of radiation and matter domination dark energy fields are slow rolling and 

(/)/--Ai/3i/ , ^2'-X2/SH , (68) 

where prime denotes derivative respects to the physical time. For the whole system in the quintessence phase (f)i will 
have a larger kinetic energy, and in the radiation dominant epoch 

H = l/2t , x' = ±lt^Jx^-Xl . (69) 

On the other hand from Ea. ((66j) we have 

vM)x' = v^A<l^iW, + v^,{cj,2W2 , (70) 

combining Eas. (|68l69l7()|) we can easily get 

VAx)=T^/)^2, (71) 
consequently the effective potential of x analytically is 

v{x) = ±^J^l-^Ux-xo) , (72) 

where xo can be easily set by the initial conditions of 0^ and the sign of or " is somewhat optional. The 
arguments above applies for the case when the total EOS of the system is restricted to be no larger than —1, the 
effective scalar will behave like phantom. 

On the evolution of perturbations we can see from Eas. (|27l28|) that the phantom and the quintessence fields obey 
the same equations. As shown in Section III although generically the two field model would have non- vanishing 
isocurvature perturbations we can choose suitable initial conditions so that the isocurvature contributions can be 
safely negligible. In this sense when the total EOS does not evolve across minus unity, the whole system can be 
equally described by an adiabatic field: both the background evolution and the adiabatic perturbations, as shown 
similarly in Refs.fT^L ItqL l8flj in the infiationary universe. 

We have demonstrated in the previous paragraphs the equivalence between the two-field quintom model and the 
single scalar field model when the EOS of the system does not cross over —1. However if the total equation of state 
for the double fields does cross over —1, this system will not be able to be described effectively by a single scalar field. 
To study the perturbations of the dark energy models with EOS across —1, we introduce a small positive constant c 
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to divide the whole region of the ahowed value of the EOS w into three parts: 1) w > — 1 + c; 2) —l-\-c>w > — 1— c; 
and 3) w < —1 — c. For 1) the EOS is always larger than —1 and for 3) w is always less than —1. For both cases the 
system with two fields as shown above can be described effectively by a single scalar field with a potential satisfying 

m 



dx^ 2^ ^ 



d ^ ,9 / 7 3 



a \2 2 



4(1 + ^) 2 



w 

- + wH{'iw + 2) 



(73) 



where is for the case 1) and " for the case 3). One can see that is divergent and there would be a 
discontinuity in the derivative V at the turning point oi w = —1, which corresponds to c ^ 0. As an example in 
Fig. 0we give the reconstructed potential of the effective x fi^ld for an oscillating quintom. One can see x behaves 
like quintessence when w > —1 and like phantom when w < —1. The reconstructed potential is well defined except 
in region 2) when the EOS gets across —1, where there is a sharp discontinuity on V'{x)- 

For the case 2), different from those in 1) and 3), the perturbations cannot be fully described by a single adiabatic 
field. However as we learn from the above, for the realistic quintom models the perturbations in the region 2) will 
be continuous and not divergent, i.e. 5 and are continuous, and the derivatives of 5 and are finite. A good 
approximation to the perturbation in region 2) is requiring it to match to the regions 1) and 3) at the boundary. 
Practically we take 5 and to be constant matching to regions 1) and 3) at the boundary and set 

j = , ^ = 0. (74) 

In the numerical calculation the constant c is a very small number, the approximation above lies in a very close 
neighborhood oi w = —1. In practice in our numerical calculations we've limited the range to be |Ak; = c| < 10~^. 
Since the region 2) is extremely limited, neglecting the evolutions of perturbations as shown in (74) is quite safe and 
well approximated. Thus we can use Eqs.(|S2 03 to study the effects of perturbations in models with parametrized 
EOS. We have also numerically checked the validity of Eq.(74) and found their contributions to the observed CMB 
and LSS power spectra are very small. The procedure of our checking is listed as follows: 

1. Start with the two-field model of quintom and record w{a)^ compute CMB and LSS spectra with perturbations. 

2. Build a code in to include dark energy perturbations with parametrized EOS. Include perturbations 
by setting Eq.JT^ and treating 5 .,9 as continuous. 

3. Interpolate w(a) in the code with parametrized EOS, compute CMB and LSS spectra and make comparisons 
with the results from step 1. 

With this procedure we have considered a model of oscillating quintom and found the difference is no more than 10~^, 
which is safely negligible. 

As examples now we study the effects of perturbations for several models with parametrized EOS in light of WMAP 
and SN data. The first example is given by Ref. Q , where w is parametrized by 

iL' (In a) = iL'a + '^0 cos [iL'i In (a/ttc)] (75) 

with a being the scale factor. This model has a nice feature of unifying the early infiation and the current accelerated 
expansions. In Ref.^42] the period of oscillation has been set as long as ~ 200 e- folds. It is interesting to study the 
consequences with a shorter period. Here for illustrations we fix k;^ = — 1, =20 and ac = 1. In the upper panels 
of FigOwe show the illustrative fittings when with and without the perturbations. We can see the parameter space 
has been enlarged a lot when including the contributions of the perturbations. For a second example we parametrize 
w as 



wihia) = Wa ^ w^aco^lwia ^ ac]- (76) 

In the numerical calculation we've fixed Wa = —1, wi = 50 and ac = 0. In the lower panels of Fig|Slwe can see the 
effects are still very prominent both in the separate and combined constraints, although not as strong as the example 
in (75). For the third example we take w to be non-oscillatory: 

w = wo/{l-\na) , (77) 

where the original form was firstly proposed in Ref. 82] . We find in our case SN constraints are very weak due to the 
fixed background parameters, the la regions have not been affected much by the perturbations, but the 2, 3cr regions 
have been enlarged significantly when the perturbations are taken into account. 



13 



VI. CONCLUSIONS 

In this paper we have studied the perturbations of the dynamical quintom model of dark energy in a self-consistent 
way. It is physically significant for the inclusion of quintom perturbations, both on the theoretical grounds of model 
buildings and on the fittings to the observations. Due to the singularities and instabilities of perturbations at the 
cosmological constant boundary, we have shown a new method regarding the impossibility of k-essence as a viable 
quintom model. In the realistic quintom model buildings one must include the perturbations. In general one needs 
to add extra degrees of freedom to realize the model of quintom. In the two-field model and the model with a 
d'Alembertian operator the isocurvature contributions may be safely negligible in the simplest case. We have con- 
sidered the implications of quintom perturbations on the observations of CMB and LSS. We have shown that the 
parameter space is different when one includes the perturbations of dark energy or not. In trying to constrain dark 
energy in a model independent way we have also proposed a method to include the perturbations for models of dark 
energy with parametrized EOS across —1. With some specific examples of the parametrized EOS, we show that 
the parameter space which characterizes the properties of the model will get enlarged in general when including the 
perturbations. This will lead to important consequences in the phenomenological studies on the cosmological imprints 
of dynamical dark energy, including the model of quintom. A thorough investigation of current constraints on the 
quintom model of dark energy where dark perturbations are taken into account is beyond the scope of current paper 
and will be presented elsewhere [23, IS4 • 

Overall, a dynamical quintom model is favored by current SN data and not ruled out by the combined observational 
constraints. There are still some inconsistencies today among different observations in the precision cosmology and 
the concordance ACDM model has not yet fitted well to the observations in a high enough confidence level, in this 
sense we might not be adopting the Ockham's razor with a cosmological constant. When we start from a ACDM 
model in the probe of our universe we cannot achieve more subtle physics beyond that. This is necessary to bear in 
mind for us to understand the nature of dark energy with the accumulation of the observational data. 
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FIG. 1: Effects of the two-field oscillating quintom on the observables. The mass of the phantom field is fixed at 2.0 x 10~^^ nipi 
and the mass of the quintessence field are 10 ~^°mpi (thicker line) and 4.0 x 10~^°mpi (thinner line) respectively. The upper right 
panel illustrates the evolution of the metric perturbations $ + ^ of the two models when with (solid lines) and without (dashed 
lines) dark energy perturbations. The scale is /c ~ 10~^ Mpc~^. The lower left panel shows the CMB effects and the lower 
right panel delineates the effects on the matter power spectrum when with (solid lines) and without (dashed lines) dark energy 
perturbations. 




FIG. 2: 3a WMAP and SN constraints on two-field quintom model shown together with the best fit values, mq and mp stand 
for the mass of quintessence and phantom respectively. We have fixed mp ^ 6.2 x 10~^^mpi and varied the value of tuq. 
Left panel: Separate WMAP and SN constraints. The green(shaded) area is WMAP constraints on models where dark energy 
perturbations have been included and the blue area(contour with solid lines) is without dark energy perturbations. Right panel: 
Combined WMAP and SN constraints on the two-field quintom model with perturbations (green/shaded region) and without 
perturbations(red region/contour with solid lines). 
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FIG. 3: Effects of the two-field quintom model where Weff — —1 compared with cosmological constant in CMB(WMAP), the 
metric perturbations $ + ^(the scale is k ^ 10~^ Mpc~^) and the linear growth factor. The binned error bars in the upper 
right panel are WMAP TT and TE data 74]. 
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FIG. 4: Reconstructions of the effective adiabatic single field x ^he framwork of oscillating quintom. The background 
parameters have been chosen as m^i = 2 x 10~^°mpz, = 10~^"^mpz, initial values are (^u = O.OQmpz, = 0A5mpi and 
4>ii = 02i = early in radiation domination epoch and for this example we have Q^i = 0.2, ^ 0.54, h ^ 0.68. The red lines 
are the total EOS of dark energy and the blue lines are the total potential. The dashed lines show the cosmological constant 
boundary. The upper panel delineates the late evolutions of the EOS and potential of dark energy and the lower panel shows 
the reconstructed values of x ctnd its potential, x is a quintessence/phantom scalar when w is above/below the dashed line. 
See the text for details. 
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FIG. 5: 3a WMAP and SN constraints on the parametrized quintom models w = —1 -\- i(;o cos(k;i In a) and w — —1 -\- 
woacos{wia), shown together with the best fit values. On the left panels the green(shaded) areas are WMAP constraints on 
models where dark energy perturbations have been included and the blue areas (contours with solid lines) are without dark 
energy perturbations. On the right panels models with perturbations are delineated in green(shaded) regions and the red 
regions (contour with solid lines) without perturbations. For illustrations we have fixed wi — 20 in the upper panels and 
It; 1 = 50 in the lower panels. 
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FIG. 6: 3(7 WMAP and SN constraints on the parametrized quintom model w — i6;o/(l — Ina) shown together with the best fit 
values. Left panel: Separate WMAP and SN constraints. The green(shaded) area is WMAP constraints on models where dark 
energy perturbations have been included and the blue area( contour with solid lines) is without dark energy perturbations. 
Right panel: Combined WMAP and SN constraints on the parametrized quintom model with perturbations (green/shaded 
region) and without perturbations(red region/contour with solid lines). 



